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We formulate a theory of spin dependent transport of an 
electronic circuit involving ferromagnetic elements with non- 
coUinear magnetizations which is based on the conservation 
of spin and charge current. The theory considerably simpli- 
fies the calculation of the transport properties of complicated 
ferromagnet-normal metal systems. We illustrate the theory 
by considering a novel three terminal device. 
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Electron transport in hybrid systems involving ferro- 
magnetic and normal metals has been shown to exhibit 
new phenomena due to the interplay between spin and 
charge. The giant magnetoresistance (GMR) effect in 
metallic magnetic multilayers is a result of spin depen- 
dent scattering ||l| . The manganese oxides exhibit a colos- 
sal magnetoresistance g] due to a ferromagnetic phase 
transition. The dependence of the current on the rela- 
tive angle between the magnetization directions has been 
reported in transport through tunnel junctions between 
ferromagnetic reservoirs P|. Transport involving ferro- 
magnets with non- coUinear magnetizations has also been 
studied theoretically in Ref. [Q 

Johnson and Silsbee demonstrated that spin dependent 
effects are also important in systems with more than two 
terminals [pi. Their ferromagnetic- normal- ferromagnetic 
(F-N-F) device manifests a transistor effect that depends 
on the relative orientation of the magnetization direc- 
tions. Recently another three terminal spin electronics 
device was realized; a ferromagnetic single-electron tran- 
sistor 1^ . In this case the current depends on the relative 
orientation of the magnetization of the source, the island 
and the drain and of the electrostatic potential of the 
island tuned by a gate voltage . 

These examples illustrate that devices with ferromag- 
netic order deserve a thorough theoretical investigation. 
Inspired by the circuit theory of Andreev reflection ||] 
we present a finite-element theory for transport in hybrid 
ferromagnetic-normal metal systems based on the conser- 
vation of charge and spin current. We demonstrate that 
spin-transport can be understood in terms of 4 general- 
ized conductances for each contact between a ferromag- 
net and a normal metal. The relations between these 
conductance parameters and the microscopic details of 
the contacts are derived and calculated for diffuse, tun- 
nel and ballistic contacts. Finally, we illustrate the the- 



ory by computing the current through a novel 3-terminal 
device. 

Let us first explain the basic idea of the finite-element 
theory of spin-transport. The system can be divided into 
(normal or ferromagnetic) "nodes" , where each node is 
characterized by the appropriate generalization of the 
distribution function, viz. a 2 x 2 distribution matrix 
in spin space. The nodes are connected to each other 
and to the reservoirs by "contacts" which limit the total 
conductance but are arbitrary otherwise. The charge and 
spin current through the contacts is related to the dis- 
tribution matrices of the adjacent nodes. Provided these 
relations are known, we can solve for the 2x2 distribution 
matrices in the nodes under the constraint of conserva- 
tion of spin and charge current in each node and thus 
determine the transport properties of the system. These 
macroscopic relations for each contact can be found in 
terms of the microscopic scattering matrices in the spirit 
of the Landauer-Biittiker formalism g| . The scattering 
matrices can be calculated using different models like a 
two-spin band model or realistic band-structures and for 
various contacts, e.g. ballistic or diffuse wires or tun- 
nel junctions. Phase coherent scattering as in a resonant 
tunneling devices and effects like the Coulomb blockade 
can be included in principle by calling the double bar- 
rier a "contact" with complex scattering properties, but 
these complications will be disregarded in the following. 

The device depicted in Fig. |^ will serve to illustrate our 
approach. Several contacts attach a normal metal node 
to (ferromagnetic or normal) metallic reservoirs. We as- 
sume that the resistances of the contacts are much larger 
than the resistance of the node. This is fulfilled when 
the area of the contact is sufficiently smaller than the 
cross-section of the node or when the contacts are in the 
tunneling regime. The current through the system and 
the distribution matrix in the node are determined by 
the properties of the contacts. The reservoirs are sup- 
posed to be large and in local equilibrium with a chemi- 
cal potential /Uq, where the subscript a labels the reser- 
voirs. The energy dependent distribution matrix in the 
(ferromagnetic) reservoir is then diagonal in spin-space 
/^(e) = l/(e,/ia), where hat (") denotes a 2 x 2 matrix 
in spin space, 1 is the unit matrix and /(e, /Iq) is the the 
Fermi-Dirac distribution function. The direction of the 
magnetization is denoted by the unit vector hIq. When 
the chemical potentials of the reservoirs are not identical, 
the normal metal node is not in equilibrium and there can 
be a spin- accumulation on the normal metal node. The 



distribution is therefore represented by a 2 x 2 matrix in 
spin-space, /^(e) which allows a spin accumulation with 
arbitrary direction of the spins. The normal metal node 
is considered to be large and chaotic either because of im- 
purity scattering inside the node or because of scattering 
at irregularities of its boundary. The distribution matrix 
inside the node is therefore isotropic in momentum space 
and depends only on the energy of the particle. 

The current through a contact is determined by its 
scattering matrix, the Fcrmi-Dirac distribution function 
of the adjacent ferromagnetic reservoir and the 2x2 non- 
equilibrium distribution matrix in the normal node. The 
current is evaluated close to the contact on the normal 
side. The 2x2 current in spin-space per energy interval 
at energy e leaving the node is 
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where M is the number of propagating channels, f"™(e) 
is the reflection matrix for an electron coming from the 
normal metal in mode m being reflected to mode n and 
i"'"(e) is the transmission matrix for an electron from the 
ferromagnet in mode m transmitted to the normal metal 
in mode n. The total current is obtained by integrat- 
ing over the energies, / = J dei{e). The current in the 
contact is thus completely determined by the scattering 
matrix of the contact, and the distribution matrices. 

The 2x2 non-equilibrium distribution matrix in the 
node in the stationary state is uniquely determined by 
current conservation 
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where a labels different contacts and the term on the 
right hand side describes spin relaxation in the normal 
node. The right hand side of Eq. (g) can be set to zero 
when the spin current in the node is conserved, i.e. when 
an electron spends much less time on the node than the 
spin-flip relaxation time Tgf. If the size of the node in 
the transport direction is smaller than the spin-flip dif- 
fusion length ^sf = v^DTsf, where D is the diffusion coef- 
ficient then the spin relaxation in the node can be intro- 
duced as (9/^/9t)rci = (iTr(/^)/2-/^)/Tst. If the size 
of the node in the transport direction is larger than Igf 
the simplest finite-element transport theory fails and we 
have to use a more complicated description with a spa- 
tially dependent spin distribution function |Q. Eq. (g) 
gives the 2x2 distribution matrix of the node in terms 
of Fermi-Diract distribution functions of the reservoirs. 
These distribution functions are determined by voltages 
of the reservoirs. Those voltages are either set by voltage 
sources or determined by conventional circuit theory. 

We will now demonstrate that the relation (|^) be- 
tween the current and the distributions has a general 



macroscopic form. Spin-flip processes in the contacts are 
disregarded, so that the reflection matrix for an incom- 
ing electron from the normal metal can be written as 
f"™ = X^s"*''"™' where s =t,i, ''"'" ^^^ the spin de- 
pendent reflection coefficients in the basis where the spin 
quantization axis is parallel to the magnetization in the 
ferromagnet, u^ — (l-|-(T-m)/2, u^ = (1 — (T-m)/2 and & 
is a vector of Pauli matrices. Similarly for the transmis- 
sion matrix i'^^(i^^)* = ^^ u^\Vl'^\^ , where i^" are the 
spin dependent transmission coefficients. Using the uni- 
tarity of the scattering matrix, we find that the general 
form of the relation P) reads 

-G^^u^f'^u^ ~ {G^^yu^f'^u^ , (3) 

where we have introduced the spin dependent conduc- 
tances G^ 
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and the mixing conductance 
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We thus see that the relation between the current 
through a contact and the distribution in the ferromag- 
netic reservoir and the normal metal node is determined 
by 4 conductances, the two real spin conductances (G^, 
G^ ) and the real and imaginary parts of the mixing con- 
ductance G^^. These contact-specific parameters can 
be obtained by microscopic theory or from experiments. 
The spin conductances G^ and G^ have been used in de- 
scriptions of spin-transport for a long time [y . The mix- 
ing conductance is a new concept which is relevant for 
transport between non-coUinear ferromagnets. The mix- 
ing conductance rotates spins around the magnetization 
axis of the ferromagnet. Note that although the mix- 
ing conductance is a complex number the 2x2 current 
in spin-space is hcrmitian and consequently the current 
and the spin-current in an arbitrary direction given by 
Eq. (H) are real numbers. Generally we can show that 
ReG^T > (G^ -|- G^)/2. Below we present explicit results 
for the conductances when the contacts are in the diffuse, 
tunneling and ballistic regimes. 

For a diffuse contact Eq. (0) can quite generally be 
found by the Green function technique developed in Ref. 
pT[ . Here we use a much simpler approach based on 
the diffusion equation. On the normal metal side of the 
contact the boundary condition to the diffusion equation 
is set by the distribution matrix in the node / . On 
the ferromagnet side of the contact the boundary con- 
dition is set by the equilibrium distribution function in 
the reservoir f^l. In a ferromagnetic metal transport 



of spins non-collinear to the local magnetization leads 
to a relaxation of the spins since electrons with differ- 
ent spins are not coherent. This causes an additional 
resistance, which as other interface related excess resis- 
tances, is assumed to be small compared to the diffuse 
bulk resistance. Sufficiently far from the ferromagnetic- 
normal metal interface the distribution function of the 
electronic states in the ferromagnet can always be repre- 
sented by two components. Only the spin-current paral- 
lel to the magnetization of the ferromagnet is conserved. 
We denote the cross-section of the contact A, the length 
of the ferromagnetic part of the contact L^ , the length 
of the normal part of the contact L^, the (spin depen- 
dent) resistitivity in the ferromagnet p^", the resistivity 
in the normal metal p^, so that the (spin dependent) 
conductance of the ferromagnetic part of the contact is 
qDFs _ A/{p^''L^) and the conductance of the normal 
part of the contact is G^^ = A/[p^L^). Solving the 
diffusion equation V^/ = on the normal and ferromag- 
netic side with the boundary conditions above, we find 
the current through a diffuse contact: 
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where the total spin dependent conductance is l/G^" = 
IjQDFs _^ i/G^iw, This resuh can be understood as 
a specific case of the generic Eq. (0) with G^ — G^\ 
G^ = G^^ and G^^ = G^^. The mixing conductance 
in the diffuse limit therefore depends on the conductance 
of the normal part of the contact only, which is a con- 
sequence of the relaxation of spins non-collinear to the 
magnetizations direction in the ferromagnet. 

For a ballistic contact, we use a simple semiclassical 
model proposed in Ref. ||l^. In this model the chan- 
nels are either completely reflected or transmitted, with 
N"^ and N-^ being the number of transmitted channels 
for different spin directions. Substituting this in (|l|) 
we find that the spin conductance G^^ = {e'^/h)N' , 
qBI __ (^^'2 /fi^]\[l and the mixed conductance is de- 
termined by the lowest number of reflected channels, 
G^Ti = max(G^^ G^^) and is real. 

For a tunneling contact we can expand Eq. (1^) in 
terms of the small transmission. We find that ReG™^ = 
(G^T + G^^)/2, where G^T and G^^ are the tunneling 
conductances. The imaginary part of G'^^^ can be shown 
to be of the same order of magnitude as G^^ and G'^^ 
but it is not universal. 

We will now illustrate the theory by computing the 
current through the 3-terminal device shown in Fig. 0. 
A normal metal node (N) is connected to 3 ferromag- 
netic reservoirs (Fl, F2 and F3) by arbitrary contacts 
parameterized by our spin-conductances. A source-drain 
bias voltage V applied between reservoir 1 and 2 causes 
an electric current / between the same reservoirs. The 



charge flow into reservoir 3 is adjusted to zero by the 
chemical potential 113. Still, the magnetization direction 
m3 influences the current between reservoir I and 2. We 
assume that spin relaxation in the normal node can be 
disregarded so that the right hand side of (||) is set to 
zero. Furthermore, we assume that the voltage bias V 
is sufficiently small so that the energy dependence of the 
transmission (reflection) coefficients can be disregarded. 
To further simplify the discussions the contacts 1 and 2 
are taken to be identical, gJ = gJ ee G^, G| = G^ = G^ 



^Ti 



and G[^ = G 

the conductances GJ, G3 and GJ 



G^K Contact 3 is characterized by 
We find the distri- 
bution in the normal node by solving the 4 linear Eqs. 
(0). The current through the contact between reservoir 

1 (2) and the node is obtained by inserting the resulting 
distribution for the normal node into Eq. (@). 

When the magnetizations in reservoir 1 and 2 are par- 
allel there is no spin-accumulation since contacts 1 and 

2 are symmetric and consequently ferromagnet 3 does 
not affect the transport properties. The current is then 
simply a result of two total conductances G = G^ -I- G-^ 
in series, / — GV/2. The influence of ferromagnet 3 is 
strongest when there is a significant spin accumulation in 
the normal metal node, and in the following the magneti- 
zations of the source and drain reservoirs are antiparallel, 
mi • m2 = — 1. We denote the relative angle between the 
magnetization in reservoir 3 and reservoir 1 (reservoir 2) 
^3 (n — 63). The current is an even function of ^3 and 
symmetric with respect to ^3 — > tt — ^3 as a result of 
the symmetry of the device, e.g. the current when the 
magnetizations in reservoir 1 and 3 are parallel equals 
the current when the magnetizations in reservoir I and 3 
are antiparallel. Due to the finite mixing conductance at 
non-collinear magnetization the third contact acts as a 
drain for the spin-accumulation in the node, thus allow- 
ing a larger charge current between reservoir 1 and 2. The 
relative increase of the current due to the reduced spin- 
accumulation A3(6l3) = [/(6I3) - /(6I3 = O)]/I{03 = 0), 
is plotted in Fig. H as a function of 63. The maximum 
of A3 is achieved at 6*3 ~ 7r/2 (6*3 ~ 3tt/2) and equals 

0) 
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introducing the total conductance of the contact Gi = 
G] + GJ, the polarization of the contact Pi = (GJ — 
GI)/{GI + GJ) and the relative mixing conductance 
7^ = 2G|V(G'I + GJ)- The influence of the direction of 
the magnetization of the reservoir 3 increases with in- 
creasing polarization P and increasing relative mixing 
conductance 773 and reaches its maximum when the total 
conductances are of the same order G3 ^ G. Note that 
the physics of this three terminal device is very different 
from that of Johnson's spin transistor J5|; the latter oper- 
ates with coUinear magnetizations of two ferromagnetic 



contacts whereas the third may be normal. 

In conclusion we have proposed a finite-element trans- 
port theory for spin transport in mesoscopic systems. In 
the presence of ferromagnetic order a contact can be de- 
scribed by 4 conductance parameters which we obtained 
explicitly for diffuse, ballistic and tunnel contacts. We 
have applied the theory to a novel three terminal device 
with arbitrary contacts. 
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FIG. 1. A normal node connected to ferromagnetic reser- 
voirs characterized by the chemical potentials /!„ and the 
magnetization vector nic, . The distribution matrix in the nor- 



mal node / 
currents ia- 



can be found from the Kirchoff rules for the spin 



FIG. 2. The 3-terminal device where a normal metal node 
is connected to ferromagnetic reservoirs. An applied bias 
causes a source-drain current between Fl and F2. The charge 
current into F3 is adjusted to vanish by /is. The magnetiza- 
tion direction of ferromagnet F3 controls the current. 



FIG. 3. The current Tr(ii) = -Tr(i2) in the three ter- 
minal device as a function of the magnetization direction 
^3. The straight, dashed and dotted lines correspond to the 
relative mixing conductance Re»72 = 2, 5, and 10 respec- 
tively. The other parameters were set to P = 0.4, P3 — 0.1, 
G — G3,Rer]i — 1.0, Imr; = = Imr/s. 
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